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This  paper  presents  an  enhanced  nonlinear  PID  (NPID)  controller  to  follow  a 
preselected  speed  profile  of  brushless  DC  motor  drive  system.  This  objective 
should  be  achieved  regardless  the  parameter  variations,  and  external 
disturbances.  The  performance  of  enhanced  NPID  controller  will  be 
investigated  by  comparing  it  with  linear  PID  control  and  fractional  order  PID 
(FOPID)  control.  These  controllers  are  tested  for  both  speed  regulation  and 
speed  tracking.  The  optimal  parameters  values  of  each  control  technique 
were  obtained  using  Genetic  Algorithm  (GA)  based  on  a  certain  cost 
function.  Results  shows  that  the  proposed  NPID  controller  has  better 
performance  among  other  techniques  (PID  and  FOPID  controller). 
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1.  INTRODUCTION 

Recently  Brushless  DC  (BLDC)  motors  have  been  replaced  conventional  DC  motors  in  many 
industrial  applications  such  as  steel  rolling  mills,  electric  trains,  electric  automotive,  aviation  and  robotics  [1]. 
This  because  BLDC  motor  has  many  advantages  compared  to  their  DC  counter -part  such  as  long  operating 
life,  high  dynamic  response,  high  efficiency  and  higher  speed  range  [2].  BLDC  motor  is  driven  by  DC 
voltage  source  but  current  commutation  is  done  by  solid-state  switches  [3].  The  commutation  instants  are 
determined  by  the  rotor  position  which  is  detected  by  position  sensors  or  by  using  sensorless  techniques  [4]. 

Proportional-Integral-Derivative  (PID)  controllers  have  been  used  for  several  decades  in  industries 
for  process  control  applications  [5].  The  reason  for  their  wide  popularity  lies  in  the  simplicity  of  design  and 
the  good  performance  including  the  low  percentage  overshoot  and  the  small  settling  time  for  slow  process 
plants  [6].  The  performance  of  PID  controllers  can  be  further  improved  by  using  of  the  fractional  order 
derivatives  and  integrals. 

The  fractional  order  PID  controllers  include  two  more  parameters  where  the  (A,)  and  (p)  are  the 
power  of  (s)  in  integral  and  derivative  actions  respectively  beside  the  three  well  known  parameters 
proportional  ( KP ),  integral  (Kj)  and  derivative  (KD)  parameters  [7].  These  additional  parameters  increase  the 
flexibility  and  robustness  of  this  controller,  and  hence  enhancing  its  dynamic  performance  compared  to  its 
integer  counterpart  [8], [9]. 

Both  the  PID  and  FOPID  controller  still  have  fixed  gains  which  are  not  enough  to  deal  with  high 
complicated  dynamic  systems  [10],[1 1].  It  is  known  that  most  of  dynamic  systems  has  nonlinear  behavior  so, 
the  linear  PID  control  is  not  adequate  to  achieve  high  performance.  As  example,  the  BLDC  motor  consists  of 
nonlinear  components  such  as  rectifier  circuit  and  six  step  inverter  [12].  This  study  presents  a  new  proposed 
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nonlinear  PID  (NPID)  controller  which  contain  on  nonlinear  gains  incorporated  with  the  fixed  gains  of  PID 
controller  [12].  These  nonlinear  gains  enjoy  the  advantage  of  high  initial  gain  to  attain  a  fast  dynamic 
response,  followed  by  a  low  gain  to  prevent  an  oscillatory  behavior  [13]. 

The  performance  of  the  PID,  FOPID  and  NPID  controller  are  highly  affected  by  selecting  its 
parameters  to  meet  the  desired  performance  of  user  for  a  given  process  plant. 

It  has  been  reported  in  several  of  researches  that  use  the  optimization  techniques  to  find  the  proper 
parameters  of  controller  is  very  useful  and  reliable  compared  to  traditional  techniques  (Ziegler-Nichols 
method).  There  are  different  types  of  optimization  techniques  such  as  Genetic  Algorithm  (GA),  particle 
swarm  and  ant  colony  [2] .  In  this  paper,  the  GA  will  be  used  to  find  the  proper  controller  parameters  of  each 
control  technique  according  to  a  certain  cost  function. 

The  main  contribution  of  paper  is  developing  a  new  technique  for  nonlinear  PID  (enhanced  NPID) 
control  and  compares  it  with  other  control  techniques  (PID  and  FOPID  control)  to  confirm  the  precedence. 
The  paper  is  organized  as  follows:  Section  II  presents  the  dynamic  modeling  of  BFDC  motor  drive  system. 
The  control  techniques  are  included  in  Section  III.  Section  IV  provides  the  simulation  results.  Finally, 
Section  V  concludes. 


2.  DYNAMIC  MODEL  OF  BLDC  MOTOR 

In  order  to  study  the  different  proposed  control  techniques  on  a  BLDC  motor,  it  is  essential  to 
introduce  the  used  dynamic  model  of  this  motor  The  transfer-function  based  mathematical  models  are  widely 
used  in  automatic  control  fields.  Some  control  design  and  analysis  methods,  such  as  the  root-locus  method 
and  the  frequency-response  method,  are  also  developed  based  on  the  system  transfer  function 

The  transfer  function  of  the  BLDC  motor  at  no  load  may  be  written  as  follows  [14], [5]: 
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Where  Ud  is  the  DC  bus  voltage,  tm  is  the  mechanical  time  constant,  te  is  the  electrical  time 
constant,  R  and  L  are  the  terminal  phase  to  phase  resistance  and  inductance,  J  is  the  rotor  moment  of  inertia, 
0)  is  the  rotor  speed,  Ke  is  the  coefficient  of  phase  back-EMF  constant,  and  KT  is  the  coefficient  of  phase 
torque  constant. 

Figure  1  shows  the  block  diagram  of  the  proposed  drive  system.  It  consists  of  a  BLDC  motor 
equipped  with  a  hall-effect  sensor  which  determines  the  rotor  position  and  sends  a  feedback  signal  to  the 
control  unit.  The  control  unit  consists  of  a  microcontroller  and  a  driving  circuit.  The  microcontroller  uses  the 
feedback  signals  of  the  motor  position  and  delivers  the  control  pulses  to  the  driving  circuit.  Then,  the  driving 
circuit  isolates  between  the  control  and  power  circuit  as  well  as  drive  the  proper  pulses  for  the  inverter  power 
circuit.  The  used  inverter  type  is  a  voltage  source  inverter.  It  converts  the  DC  supply  to  a  suitable  one  in 
order  to  drive  the  BLDC.  The  DC  supply  is  achieved  by  using  an  uncontrolled  rectifier  followed  by  a  DC  link 
capacitor  [15]. 


Figure  1.  Brushless  DC  motor  drive  system 
The  parameters  of  the  BLDC  motor  are  listed  in  Table  1. 
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Table  1.  The  BLDC  Motor  Parameters. 


Rating 

Symbol 

Value 

Units 

DC  resistance 

R 

0.57 

O 

Inductance 

L 

1.5 

mH 

Torque  constant 

Kt 

0.082 

N.m/A 

No.  of  Poles 

P 

4 

Rated  torque 

T 

ip 

0.42 

N.m 

Rated  Voltage 

V 

36 

V 

Rotor  Inertia 

J 

23e-6 

Kg.  m2 

Friction  coefficient 

Bv 

0.0000735 

N.M.S 

Rated  Speed 

O) 

4000 

RPM 

Rated  current 

I 

5 

A 

Figure  2  demonstrates  the  open  loop  response  of  BLDC  motor  drive  system  model.  It  can  be  noted 
that  at  time  0.1  sec  the  BLDC  motor  exposed  to  sudden  load  correspond  50%  of  rated  torque,  the  rotor  speed 
will  be  reduced  to  3200  RPM.  Figure  3  shows  the  corresponding  phase  current  of  BLDC  motor.  It  can  clear 
that  the  phase  current  will  be  increased  at  time  0.1  sec  to  ±2.5  A. 


Figure  2.  Open  loop  response  of  BLDC  motor 
drive  system  model 


Figure  3.  The  corresponding  phase  current  of  BLDC 
motor  model 


3.  CONTROL  TECHNIQUES 

This  section  illustrates  the  theory  of  three  different  control  techniques.  The  first  technique  is  the 
conventional  PID  control  while,  the  second  technique  is  the  fractional  order  PID  control.  The  third  technique 
is  the  nonlinear  PID  control.  The  parameters  of  each  control  technique  were  obtained  by  GA  optimization 
technique  based  on  a  certain  cost  function. 

3.1.  PID  Control 

Proportional-integral-derivative  (PID)  controllers  have  been  the  most  famous  and  the  most 
commonly  used  industrial  controllers  in  the  previous  years.  The  reputation  and  widespread  use  of  PID 
controllers  are  attributed  primarily  to  their  simplicity  and  performance  characteristics  [16].  PID  controllers 
have  been  applied  for  control  of  various  dynamical  systems  ranged  from  industrial  process  to  complicated 
dynamics  such  as  ship  and  aircraft  [17]. 

The  transfer  function  of  the  PID  controller  is  K(s)  =  KP  +  —  +  Kds.  Where  Kp,  Kh  kd  are 

proportional,  integral  and  differential  gains  respectively.  The  task  of  each  part  of  a  PID  controller  can  be 
described  as  follows,  the  proportional  part  decreases  the  error  responses  of  the  system  to  disturbances,  the 
integral  part  rejects  the  steady-state  error,  and  finally  the  derivative  part  dampens  the  dynamic  response  and 
improves  the  system  stability.  The  problem  in  the  PID  controller  is  to  select  the  proper  three  parameters  to  be 
suitable  for  the  controlled  plant.  There  are  different  methods  to  find  the  parameters  of  PID  controller  such  as 
try  and  error  and  Ziegler-Nichols  method  but,  most  of  these  methods  are  rough  roads.  In  this  paper  the  GA 
optimization  technique  will  be  used  to  obtain  the  optimal  values  of  PID  controller  parameters  based  on  the 
cost  function  as  shown  in  (4). 
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This  cost  function  can  satisfy  the  designer  demand  using  the  weighting  factor  value  (P).  The  factor 
is  set  larger  than  0.7  to  decrease  over  shoot  and  steady-state  error.  If  this  factor  is  set  smaller  than  0.7  the  rise 
time  and  settling  time  will  be  decreased  [18]. 

Matlab  genetic  algorithm  toolbox  is  used  to  tune  the  PID  controller  parameters  according  to  the 
previous  objective  function.  The  first  and  the  most  critical  step  is  to  convert  the  problem  into  suitable  GA 
chromosomes  and  then  construct  the  population.  Some  works  recommend  20  to  100  chromosomes  in  one 
population.  The  higher  the  chromosomes  number,  the  better  chance  to  get  the  optimal  results.  However, 
using  more  chromosomes  requires  more  execution  time.  In  this  paper,  it  is  desired  to  use  80  chromosomes  in 
each  generation. 

Each  chromosome  comprises  of  three  parameters,  i.e.  KP)  KIf  KD  with  value  bounds  varied  depend 
on  the  used  objective  function.  In  order  to  reach  to  accurate  values  for  the  PID  parameters,  the  initial  values 
of  KP,  KIf  Kd  are  determined  using  Ziegler-Nichols  rule.  The  population  in  each  generation  is  represented  by 
80  x  4  population  (P)  matrix  as  given  by  (3),  depends  on  the  chromosomes  number  in  population,  where  each 
row  represents  a  one  chromosome  that  comprise  KP,  Kit  Kd  values  and  the  last  column  is  added  to 
accommodate  fitness  values  (F)  of  corresponding  chromosomes. 


\KV1 

Kn 

Kd  i 

Fll 

Kp2 

Ku 

Kd2 

f2 

Pl  = 

&  :  : 

Kin 

Kdn 

Fn_ 

Where  n  is  the  number  of  chromosomes. 


3.2.  FOPID  Control 

The  most  common  form  of  a  fractional  order  PID  controller  is  the  PIAD^,  including  an  integrator  of 
order  X  and  a  differentiator  of  order  p  as  demonstrated  in  Figure  4  where  p  and  X  can  be  any  real  numbers. 
The  controller  transfer  function  has  the  form  [19]: 

Gc  =J^  =  kp  +  ki^I  +  kds'i,(A,fi>  0)  (6) 


Figure  4.  FOPID  controller  block  diagram 


Figure  5.  Expanding  from  Point  to  Plane 


Conventional  PID  controller  can  be  obtained  by  adjusting  A  =  1  and  fi  —  1.  So,  the  conventional 
PID  controller  is  special  case  from  the  FOPID  controller  which  enlarges  the  PID  control  from  point  to  plane 
as  shown  in  Figure  5.  This  extension  provides  us  much  more  flexibility  and  reliability  and  provides  a  chance 
to  enhance  the  performance  of  control  system  and  that  is  through  the  two  more  parameters  of  FOPID. 
Moreover,  this  improves  the  robustness  of  the  system  against  any  unexpected  variables  and  makes  it  more 
stable  where  by  adjusting  the  two  more  parameters  of  FOPID  can  compensate  any  disturbances  in  system  . 
But,  now  the  task  of  the  designer  selects  the  five  parameters  instead  of  three  parameters  which  will  be  more 
challenging  for  achieving  a  certain  performance  [19]. 
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In  this  paper,  the  FOPID  control  parameters  will  be  adjusted  using  GA  optimization  based  on  the 
previous  cost  function.  The  first  step  is  to  encrypt  the  problem  into  appropriate  GA  chromosomes.  Each 
chromosome  includes  five  parameters,  i.e.  KP,KIf  KD ,  X  and  p  with  value  bounds  varied  depend  on  the  used 
cost  function.  The  initial  values  of  KP,  KIt  KD,  X  and  p  are  obtained  using  Ziegler-Nichols  rule  to  reach 
accurate  values.  The  population  in  each  generation  is  represented  by  80  x  6  population  (P2)  matrix  as  given 
by  (4),  depends  on  the  chromosomes  number  in  population,  where  each  row  represents  a  one  chromosome 
that  comprise  KP)  Kit  Kd,X  and  p  values  and  the  last  column  is  added  to  accommodate  fitness  values  (F)  of 
corresponding  chromosomes. 

Kv  i  Ktl  Kdl  X1  px  F1 ' 

Kp2  Ki2  Kd2  ^2  P-2  ^2 

P2  = .  (7) 

_Kpn  Kin  Kdn  Xn  pn  Fu_ 

Where  n  is  the  number  of  chromosomes. 


3.3.  NPID  Control 

For  controlling  a  nominal  physical  process,  the  demands  for  high-performance  control  with  different 
operating  conditions  or  environmental  parameters  are  often  beyond  the  capabilities  of  simple  PID  controllers. 
In  order  to  improve  the  performance  of  linear  PID  controllers,  many  techniques  have  been  established  to 
increase  the  adaptability  and  robustness  by  adopting  the  self-tuning  method,  general  predictive  control, 
neural  networks  strategy  and  fuzzy  logic,  and  other  methods  [13]. 

Among  these  techniques,  nonlinear  PID  (NPID)  control  is  presented  as  one  of  the  most  effective  and 
simple  method  for  industrial  applications.  The  nonlinear  PID  (NPID)  control  has  found  in  form  two  broad 
categories  of  applications  [10].  The  first  category  particular  to  nonlinear  systems,  where  NPID  control  is  used 
to  accommodate  the  nonlinearity  while,  the  second  category  deals  with  linear  systems,  where  NPID  control  is 
used  to  obtain  performance  not  realizable  by  a  linear  PID  control,  such  as  increased  damping,  reduced  rise 
time  for  step  or  rapid  inputs,  improved  tracking  accuracy,  and  friction  compensation. 

The  NPID  controllers  have  the  advantage  of  high  initial  gain  to  achieve  a  fast  response,  followed  by 
a  low  gain  to  avoid  an  oscillatory  behavior.  In  this  study,  a  performance  enhancement  to  the  conventional 
linear  PID  controller  is  proposed  by  combining  a  sector-bounded  nonlinear  gain  into  a  linear  fixed  gain  PID 
control  architecture  [12]. 

In  this  paper,  the  proposed  enhanced  nonlinear  PID  (NPID)  controller  consists  of  two  parts.  The  first 
part  is  a  sector  bounded  nonlinear  gain  Kn(e) while,  the  second  part  is  a  linear  fixed-gain  PID  controller  (Kp, 
Ki  and  Kd).  The  nonlinear  gain  Kn(e)  is  a  sector-bounded  function  of  the  error  e(t).  the  previous  researches 
have  been  made  where  the  nonlinear  gain  Kn(e)  has  a  one  scalar  value. 

In  this  paper,  The  Kn(e )  acts  as  a  row  vector  can  expressed  as  Kn(e )  =  \Knl(e)  Kn 2(e)  Kn 3(e)] 
which  will  be  enhanced  the  NPID  controller. 

The  proposed  form  of  NPID  control  can  be  described  as  follows. 

U(t)  =  Kp[Knl(e).e(t)]  +  Kt  f*[Kn2(.e).e(t)]  dt  +  Kd  [ffn3(e).^]  (8) 

Where  Knl(e),Kn 2(e)  and  Kn3(e )  are  nonlinear  gains.  The  nonlinear  gains  represent  any  general 
nonlinear  function  of  the  error  e  which  is  bounded  in  the  sector  0  <  Kn(e)<  Kn(e) max.  There  is  a  broad 
range  of  options  available  for  the  nonlinear  gain  Kn(e ).  One  simple  form  of  the  nonlinear  gain  function  can 
be  described  as. 


Kni(e)  =  ch(wte ) 


exp(wje)+exp  (-w^e) 
2 


(9) 


Where  i  =  1, 2,  3. 


e  = 


e 


1^1  —  emax  | 

\e\  >  emax  ) 


(10) 
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Where  wt  and  emax  are  user-defined  positive  constants.  The  nonlinear  gain  Kn(e )  is  lower  bounded 
by  Kn(e) min  =  1  when  e  =  0,  and  upper-bounded  by  Kn(e) max  =  ch(w;  emax).  Therefore,  emax  stand  for  the 
range  of  variation,  and  wt  describes  the  rate  of  variation  of  Kn(e ). 

Matlab  genetic  algorithm  toolbox  is  used  to  adjust  the  NPID  controller  parameters  according  to  the 
earlier  cost  function.  The  population  in  each  generation  is  represented  by  80  x  7  population  (P3)  matrix  as 
given  by  (8).  Each  row  represents  a  one  chromosome  that  include  KP,  Kif  Kd,wlf  w2  and  w3  values  and  the 
last  column  is  added  to  adapt  fitness  values  (F)  of  corresponding  chromosomes. 


Kvi 

Kii 

Kd  i 

Wn 

W2i 

W31 

Fll 

KP2 

Ki2 

Kd2 

W12 

W22 

W32 

f2 

Kpn 

Kin 

Kdn 

Win 

Wm 

w3n 

Fn_ 

Where  n  is  the  number  of  chromosomes. 


4.  SIMULATION  RESULTS 

This  section  presents  investigates  the  performance  of  each  control  technique  where  several  of 
investigations  will  be  applied  on  the  BLDC  motor  such  as  speed  regulation,  sudden  load,  speed  tracking  and 
continuous  variable  load. 

Figure  6  demonstrates  the  speed  response  of  Both  the  NPID  control  with  scalar  and  vector  nonlinear 
gain.  It  is  obvious  that  the  rotor  of  BLDC  motor  was  exposed  to  sudden  load  by  half  rated  torque  at  time  0.1 
sec.  Also,  it  can  be  noted  that  the  performance  of  the  NPID  controller  with  vector  nonlinear  gain  has  faster 
response  than  other  kind  of  NPID  control.  Moreover,  it  can  absorb  the  disorder  rapidly  at  time  0.1  sec. 

The  rest  of  tests  will  be  performed  by  using  the  NPID  control  with  a  vector  value  and  will  be 
compared  with  GA-PID  control  and  FOPID  control. 


Figure  6.  Speed  response  of  control  techniques  during 
the  speed  regulation 


Figure  7.  Speed  response  of  control  techniques 
during  the  speed  regulation 


Figure  8  demonstrates  the  corresponding  DC  supply  current  of  each  control  technique.  It  can  be 
noted  that  the  NPID  controller  has  high  starting  current  compared  to  other  techniques.  In  case  of  sudden  load 
at  time  0.1  sec,  the  DC  supply  current  of  NPID  control  will  be  raised  rapidly  among  other  techniques  to  can 
accommodate  the  disturbance  in  a  short  time. 
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Figure  8.  The  corresponding  DC  supply  of  each 
control  technique  during  the  speed  regulation 


Figure  9.  Sinusoidal  load  torque  that  exposed  to  the 
BLDC  motor  through  the  simulation 


Table  2  summarizes  and  give  us  accurate  comparison  of  each  control  technique  performance  Figure 
6  and  7.  It  is  obvious  that  the  dynamic  response  of  GA-PID  control  is  very  slow  compared  to  other 
techniques.  Also,  the  rise  time  of  FOPID  control  close  to  the  rise  time  of  NPID  control  (with  vector  nonlinear 
gain)  but,  it  suffers  from  high  overshoot.  Moreover,  the  NPID  control  (with  vector  nonlinear  gain)  has  the 
minimum  rise  and  settling  time  and  no  overshoot. 


Table  2.  Control  Techniques  Performance 


Controller  Method 

Rise  Time  (sec) 

Settling  Time  (sec) 

Max.  Overshoot  % 

GA-PID  Control 

0.0219 

0.0342 

1.1107 

FOPID  Control 

0.0071 

0.0756 

18.13 

NPID  control  (with  scalar  nonlinear  gain) 

0.0195 

0.1301 

0.9382 

NPID  Control  (with  vector  nonlinear  gain) 

0.0067 

0.0099 

0.0001 

Sometimes,  the  BLDC  motor  exposed  to  sinusoidal  load  when  used  as  actuator  in  robotics 
applications  because  of  the  change  of  links  inertia  with  robot  movement.  So,  the  BLDC  motor  needs  robust 
control  to  can  regulate  the  speed  of  BLDC  motor  through  the  continuous  load  change  which  can  expressed  as 
sinusoidal  function  has  a  peak  0.2  N.m  as  shown  in  Figure  9. 

Figure  10  illustrates  the  speed  response  of  each  control  technique  through  the  continuous  load 
change.  It  is  noted  that  the  rotor  speed  of  BLDC  motor  oscillated  about  the  reference  speed.  Both  the  GA- 
PID  control  and  FOPID  control  have  high  deviation  about  reference  speed  while,  the  enhanced  NPID  control 
can  eliminate  the  continuous  disturbance  without  high  deviation  about  the  reference  speed. 
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Figure  10.  Performance  of  control  techniques  through 
sinusoidal  load 


Figure  1 1 .  The  corresponding  control  techniques 
output  through  sinusoidal  load 
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Table  3  demonstrates  the  deviation  values  about  the  reference  speed  of  each  control  technique.  It  is 
noted  that  the  GA-PID  control  and  FOPID  control  have  high  deviation  while,  the  NPID  control  has  very 
small  deviation  about  the  reference  speed. 


Table  3.  Control  Techniques  Deviation. 


Controller  Method 

Deviation  (RPM) 

Deviation  (%) 

GA-PID  Control 

±220 

7.34% 

FOPID  Control 

±75 

2.5% 

NPID  Control 

±5 

0.17% 

The  reference  speed  of  BLDC  motor  in  some  industrial  applications  is  not  constant  and  need  to 
change  continuously  such  as  robotics  and  electric  automotive.  So,  Figure  12  shows  the  response  of  each 
control  technique  at  different  commands  of  reference  speed  with  constant  load  correspond  10%  of  rated 
torque.  We  note  that  NPID  control  can  track  the  reference  speed  faster  than  other  techniques. 

Figure  13  illustrates  the  corresponding  DC  supply  current  of  each  controller  technique.  It  is  clear 
that  the  NPID  control  has  high  starting  current  with  each  change  in  reference  speed  compared  to  other 
techniques. 


Figure  12.  Speed  response  of  control  techniques  at 
different  commands  of  reference  speed 


Figure  13.  The  corresponding  DC  supply  current  of 
control  techniques  at  different  commands  of 
reference  speed 


5.  CONCLUSION 

In  this  paper,  an  enhanced  nonlinear  PID  (NPID)  controller  is  proposed  to  improve  the  performance 
of  BLDC  motor  and  track  a  preselected  speed/position  profile  regardless  the  parameter  variations  and  sudden 
load.  To  ensure  the  robustness  of  NPID  control  two  different  techniques  were  compared  with  it.  The  first 
technique  is  the  conventional  PID  control  while  the  second  technique  is  the  fractional  order  PID  control.  The 
genetic  algorithm  optimization  technique  is  used  to  obtain  the  proper  parameters  of  each  control  technique. 
Simulation  results  illustrate  that  the  proposed  NPID  controller  can  accommodate  the  different  types  of 
disturbance  rapidly  and  has  faster  response  compared  to  other  control  techniques. 
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